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We consider a Hamiltonian quantum theory of stationary spacetimes containing a Kerr-Newman 
black hole. The physical phase space of such spacetimes is just six-dimensional, and it is spanned 
by the mass M, the electric charge Q and angular momentum J of the hole, together with the 
corresponding canonical momenta. In this six-dimensional phase space we perform a canonical 
transformation such that the resulting configuration variables describe the dynamical properties of 
Kerr-Newman black holes in a natural manner. The classical Hamiltonian written in terms of these 
variables and their conjugate momenta is replaced by the corresponding self-adjoint Hamiltonian 
operator and an eigenvalue equation for the Arnowitt-Deser-Misner (ADM) mass of the hole, from 
the point of view of a distant observer at rest, is obtained. In a certain very restricted sense, 
this eigenvalue equation may be viewed as a sort of "Schrodinger equation of black holes". Our 
"Schrodinger equation" implies that the ADM mass, electric charge and angular momentum spectra 
of black holes are discrete, and the mass spectrum is bounded from below. Moreover, the spectrum 
of the quantity M 2 — Q 2 — a 2 , where a is the angular momentum per unit mass of the hole, is strictly 
positive when an appropriate self-adjoint extension is chosen. The WKB analysis yields the result 
that the large eigenvalues of M, Q and a are of the form \/2ri, where n is an integer. It turns out 
that this result is closely related to Bekenstein's proposal on the discrete horizon area spectrum of 
black holes. 
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I. INTRODUCTION 



Black holes are among the simplest and most beautiful objects in the universe. They are made of spacetime and 
electromagnetic field only, and they have just three classical degrees of freedom which may be taken to be the mass 
M, electric charge Q and the angular momentum J of the hole. 

Although the number of classical degrees of freedom of black holes is just three, however, one expects that there 
is an enormous number of quantum mechanical degrees of freedom associated with black holes. During some recent 
years, string theory and loop theoretic approaches to quantum gravity have greatly improved our understanding on 
the nature of these quantum mechanical degrees of freedom jjlpj . 

As it happens, there is a resemblance between black holes and hydrogen atoms. Like a black hole, a hydrogen 
atom has just three classical degrees of freedom. Indeed, the system looks very simple: An electron whirls around 
the proton, and the classical degrees of freedom may be taken to be the x—, y— and z~ coordinates of the electron. 
Quantum field theoretical investigations reveal, however, an enormous number of quantum mechanical degrees of 
freedom associated with virtual electron-positron pairs and photons. Still, the quantum mechanical properties of the 
hydrogen atom may be described, as an excellent approximation, by its non-relativistic Schrodinger equation which 
takes into account the three classical degrees of freedom only. 

The resemblance between black holes and hydrogen atoms gives rise to an interesting question of a possibility to 
construct a quantum-mechanical model of a black hole which, although it takes into account the three classical degrees 
of freedon of black holes only, nevertheless describes their quantum mechanical properties with a reasonable accuracy. 
In this paper we shall consider one such model of black holes. Of course, even classical black hole spacetimes may 
perform all sorts of vibrations and oscillations which provide them with an enormous number of additional degrees 
of freedom, but in this paper we are interested in stationary black holes only. In other words, we are quantizing 
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the stationary sector of black hole spacetimes, and such a sector is characterized by just three classical degrees of 
freedom jjj. 

Our model is based on an observation that even stationary black hole spacetimes have dynamics. More precisely, 
even stationary black hole spacetimes have a region which does not admit a timelike Killing vector field. This means 
that in a certain spacetime region the black hole spacetime geometry evolves in time no matter how we choose the 
time coordinate. It is this time evolution of black hole spacetime geometry on which we focus our attention and 
which, in our model, is responsible for the quantum mechanical properties of black holes. 

To see what this means consider, as an example, the simplest possible black hole, the Schwarzschild black hole. It 
has the spacetime metric 



ds 2 = - ( 1 - dt 2 + + t 2 (de 2 + sin 2 Odtf) . (1.1) 



One observes that when r < 2M, the coordinate r becomes timelike, and because spacetime geometry inside the 
event horizon depends on r, it evolves in time. In that region r describes the radius of the wormhole throat of the 
black hole. In a more precise manner the fact that spacetime inside the event horizon really has dynamics in its 
geometry can be seen if one considers the conformal diagram of Kruskal spacetime: When r < 2M one cannot move 
in any timelike direction without changing r, and therefore the geometry of the spacelike hypersurfaces of spacetime. 
In Reissner-Nordstrom and Kerr-Newman black hole spacetimes the dynamical region lies in the intermediate region 
between the outer and the inner horizons of the hole. 

In this paper we consider the Hamiltonian quantization of Kerr-Newman black hole spacetimes in such a manner 
that, in the classical level, the phase space coordinates of the theory describe the dynamics of the intermediate region 
between the horizons in a natural way. The Kerr-Newman solution is a specific solution to Einstein's and Maxwell's 
equations in vacuum. Because of that, we begin our investigations in Section || by a general study of the Hamiltonian 
formulation of the Einstein-Maxwell theory, paying particular attention to the boundary terms which are needed in 
asymptotically flat electrovacuum spacetimes, such as Kerr-Newman spacetimes, to make the Hamiltonian formulation 
consistent. We shall see in later sections that these boundary terms play a most fundamental role in the quantum 



theory of Kerr-Newman black holes. In Section III we calculate these boundary terms for maximally extended Kerr- 
Newman spacetimes. It turns out that from the boundary terms one can read off the mass, electric charge and angular 
momentum of the black hole. 



The study of the classical Hamiltonian dynamics of Kerr-Newman black hole spacetimes is performed in Section IV . 
Basically, our study is based on an important theorem proved by Regge and Teitelboim |Q. This theorem states, 
essentially, that the physical Hamiltonian of asymptotically flat spacetime with matter fields can be gained if we first 
solve the classical constraints, and then substitute the solutions to the constraints, in terms of the physical phase 
space coordinates of the theory, to the boundary terms at asymptotic spacelike infinity. At the first stage we take 
the phase space coordinates to be the mass M, the electric charge Q and the angular momentum J of the hole, 
together with the corresponding canonical momenta pm, Pq and pj, and we write the sum of the boundary terms, 
and hence the classical Hamiltonian, in terms of these phase space coordinates. It is unclear whether the assumptions 
of the Regge's and Teitelboim's theorem are valid for Kerr-Newman spacetime and the variables M, Q, J, pm, Pq 
and pj, but we accept this as an unproved hypothesis and see where it takes us. At the second stage we perform a 
canonical transformation from the variables M, Q and J and their canonical momenta to the new variables and their 
canonically conjugate momenta which describe better the dynamics of the intermediate region between the horizons 
of the Kerr-Newman black hole. In terms of these phase space variables we write the classical Hamiltonian of Kerr- 
Newman spacetimes in a specific foliation where the flat Minkowski time coordinate of an asymptotic observer at rest 
at a faraway infinity coincides with the the proper time of a freely falling observer at the throat of the black hole. An 
explicit example of such a foliation is presented in Appendix [b|. 

In Section ^ we proceed to quantization. A straightforward replacement of the classical Hamiltonian by the 
corresponding self-adjoint Hamiltonian operator yields an equation which, in a certain very restricted sense, may be 
considered as a sort of "Schrodinger equation of black holes". That equation is the main result of this paper. In the 
natural units, where h = c = G = 47reo = 1, and when a particular operator ordering has been chosen, it can be 
written, in terms of the configuration variables R, u and v of the theory, as: 

where VP = ^(R, it, v) is the wave function of the Kerr-Newman black hole. 

Besides of being an equation of considerable simplicity, our "Schrodinger equation" has many interesting conse- 
quences. For instance, it predicts that the mass, electric charge and angular momentum spectra of black holes are 
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discrete. In particular, the mass spectrum is bounded from below and can be made positive by means of an appro- 
priate choice of a self-adjoint extension. As a matter of fact, one may prove even more than that: It is possible to 
choose operator orderings and self-adjoint extensions such that the spectrum of the quantity 



M 2 - Q 2 - a 2 , 

where a := -j^ is the angular momentum per unit mass, is strictly positive. Regarding Hawking radiation, this is very 
interesting result: It implies that a non-extreme black hole can never become, by means of Hawking radiation, an 
extreme hole. This result is in agreement with the third law of black hole thermodynamics, and is therefore a strong 
argument in favor of the physical validity of our model. 

At the high end of the spectrum, we find that the eigenvalues of the sum of the areas of the horizons of the hole 
which we shall call, for the sake of convenience, the total area of the black hole, are of the form 

= n ■ 1671-Zf,! , (1.3) 

where n — 1, 2, 3, . . . and Zpi := J is the Planck length. As such our result is closely related, although not quite 

identical to Bekenstein's proposal ]5|-|25[| . According to that proposal, the spectrum of the outer horizon of the hole 
is of the form 

A n = n ■ 7^p! , (1.4) 

where n = 1, 2, 3, . . . and 7 is pure number of order one. Arguments in favor of the claim that it is perhaps not the 
area of the exterior horizon but the sum of the areas which should have an equal spacing in its spectrum will be given 
111 Section [v| 

Finally, our model implies a very interesting discrete spectrum for the angular momentum of Kerr-Newman black 
holes. For uncharged holes near extremality, the angular momentum eigenvalues are of the form 

J m = mh , (1.5) 

where m = 0, ±2, ±4, . . .. 

II. HAMILTONIAN FORMULATION OF EINSTEIN-MAXWELL THEORY 

The Einstein-Maxwell theory is a theory of electromagnetic field interacting with gravitational field. In this section 
we shall develope the Hamiltonian formulation of such a theory in all details, paying particular attention to the 
boundary terms appearing in asymptotically flat spacetimes as a consequence of the requirement of internal consistency 
of the theory. 

The action of the Einstein-Maxwell theory can be written, in general, as 

S = j^J d i Xy/=g~ ( ( % - F^F^ + (boundary terms) . (2.1) 

In this equation the integration is performed over the whole four-dimensional spacetime. g is the determinant of the 
spacetime metric <7 M „, and 

F/J.V := d^Ay - d v A a (2.2) 

is the electromagnetic field tensor. A^ is the electromagnetic vector potential. ^>R is the four-dimensional scalar 
curvature. 



As it is well known, we can write the action (2.1) as 

S = S| rav + + SfT + S c 9 % , (2.3) 

where 

Sl rav := J dt Jj 3 x^N(K ab K ab - K 2 + K) , (2.4a) 

S ™ := l(b / * J^^NF^ , (2.4b) 
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and 5|^ v and Sq£ are boundary terms associated with spacelike asymptotic infinities of asymptotically flat spacetimes. 
In Eqs. (2.4) the spatial integration is performed over the whole spacelike hypersurface £ of spacetime where the time 
t is constant. K ab is the exterior curvature tensor on that hypersurface, K its trace, and TZ is the three-dimensional 
scalar curvature on that hypersurface. N is the lapse function and q is the determinant of the metric q ab on the 
hypersurface S. 

The properties of the action S^ av are well known. Consider now the action Sff 1 of Eq. ( 2.4tj ). To begin with, 
consider first the case where the spacetime metric can be written as 



ds 2 



-di z + q a bdx a dx b 



(2.5) 



In other words, we have chosen time orthogonal coordinates, where the lapse N = 1, and the shift iV a vanishes 
identically. In these coordinates we can write 



SfP = I dt I d 3 x£ em 



where 



A a A b - 2A a (d b A ) + (d a A )(d b A ) 



(3) F .(% 



is the electromagnetic Lagrangian in curved spacetime. The dot means time derivative and we have defined 

{3) F ab := d a A b - d b A a , 

(3)pab ._ qaniybn l 3 )j? mn 

The canonical momentum conjugate to A a is 



P 



This relation can be inverted, and we have 



where we have defined 



47T 

A b = —p b + d b A 



Pb ■= q a bP ■ 

In terms of p a we can write the electromagnetic Lagrangian as 

'2tt 



C 



cm a 



P a A a 



q ab p a p b +P a (d a A ) 



16tt 



Hence, we get 



where 



dt I d 3 x 

s 



P a A a - n cm + Md aP a ) 



n ^^^ qabP yb + ^ )Fa p F ab 



In Eq. (2.13) we have dropped the term | J dt J E d 3 xd a (Aop a ) which can be inverted to a boundary term. 
We now include the lapse and the shift to our formulation. To include the lapse we replace dt by 

dt 1 = Ndt , 



(2.6) 



(2.7) 



(2.8a) 
(2.8b) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



and, because Aq transforms to 
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we find that for general lapse but vanishing shift the electromagnetic action is 

sgn = I dt L d3x [ paAa ~ Nn ° m + Mdapa) 

Inclusion of a non- vanishing shift is a bit more tricky. We replace dx a by 

dx' a = dx a + N a dt , 

from which it follows that Aq is replaced by 

A = Aq - N s A s . 
Moreover, at the hypersurface where x° = t + dt, A a is replaced by 



A' 



dx s 



■A s (t + dt,x b - N b dt) =A a + A a dt - (d s A a )N s dt - {d a N s )A s dt 



a dx' a ' 

Hence, we find that A a must be replaced by 

A' a = A a - (d s A a )N s - (d a N s )A s . 



(2.16) 

(2.17) 

(2.18) 
(2.19) 

(2.20) 

(2.21) 



Substituting Eqs. ( 2.19 ) and ( 2.21 ) into Eq. ( 2.17 ) we obtain an expression for the electromagnetic action in the 
presence of non- vanishing shift: 



= I dt I d"x 



where we have defined 



p a A a - NH em - N s H c s m + A (d aP a ) 



H em ._ p a (3) Fgi 



(2.22) 



(2.23) 



and we have ignored the term J dt / s d 3 xd a (A s N s p a ). 

We are now prepared to write down the whole Einstein-Maxwell action without boundary terms. The gravitational 
part S , |T av is a mere ADM action 



crgra 



where 



n sia 

Hf a 



dt / d 3 x (p ab q ab - NH siav - N s Hf av ) 



l(16ir)G abcd p ab p cd + zr^s/qll 

I 1D7T 



and 



P ab :=- — ^{K ab -q ab K) 



is the canonical momentum conjugate to q ab . 
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Gabcd '■— -p {qabqcd ~ <7ac9bd ~ 9ad9bc 



(2.24) 

(2.25a) 
(2.25b) 

(2.26) 
(2.27) 



is the Wheeler-De Witt metric. Putting the actions (2.22) and (2.24) together we get the Einstein-Maxwell action 

(2.28) 



S| rav = dt d 3 x (p ab q ab +P a A a - NH — N S H S - AqQ) , 



where 
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is the Hamiltonian constraint, 



is the diffeomorphism constraint, and 



U := H glav + H cm (2.29) 
H s := Kf av + Hf 11 (2.30) 
Q := -d aP a (2.31) 



is the Gaussian constraint. 

We shall now consider asymptotically flat spacetimes. In those kind of spacetimes we must include certain boundary 
terms, since we cannot assume the variations of the dynamical variables and their canonical momenta to vanish at 
asymptotic infinity. In what follows, we shall take the asymptotic coordinates at spatial infinity to be Cartesian 
coordinates. 

First of all, of course, we have the ADM boundary term p|j26| 

Sa£ M = - J dtN+(t)E ADM (t) , (2.32) 

where 

N+(t) := lim N(t,x a ) (2.33) 

r— >oc 

is the lapse function at the asymptotic spatial infinity, and 



is the ADM energy of spacetime. In Eq. (2.34) h mn denotes a spatial component of the linearized gravitational field 
in asymptotic Cartesian coordinates. More precisely, we have assumed spatial coordinates to become Cartesian in 
asymptotic spacelike infinity, and in these coordinates we have written the spacetime metric as g^ v = r)u V + hu V , where 
= diag(— 1, 1, 1, 1) is the flat spacetime metric. 
In addition to the ADM boundary term, which is a term associated with time evolution at asymptotic infinity, we 
have, for the non-vanish ing s hift at spatial infinity, boundary terms associated with asymptotic spatial translations. 
Variation of the action ( 2.28 ) with respect to the momentum p ab conjugate to q a b brings along a term Q 



dt J^d 3 x(N a Sp ab ) lb , 
which must be cancelled at infinity. Hence we need a boundary term 

SlT S ■=- J dtN+(t)P^ M (t) , (2.35) 

where 



is the shift at infinity, and 



N+(t) := lim N a {t,x a ) (2.36) 



-Padm := - lim 2 4 p ab dS b (2.37) 



is the ADM momentum of spacetime. 

So far we have considered terms related to pure gravity. We still have to include boundary terms related to 
electromagnetism. First of all, we observe that variation of the action with respect to the momentum p a conjugate 
to A a brings along a term 



dt I d 3 xd a (A 5p a ) 
s 



G 



which must be cancelled at infinity. Hence, we need an electromagnetic boundary term 

S%:=- J dtA+{t)Q{t) , (2.38) 

where 

At(t) := lim A {t,x a ) (2.39) 

r — >oo 

is the electric potential at infinity, and 

Q := - lim <£ Pa dS a (2.40) 



is the electric charge of spacetime. 

We are now prepared to write the whole Einstein- Maxwell action, with appropriate boundary terms. We get 



Sj: = j dt J d 3 x (p ab q ab + p a A a — NH — N S U S - A Q g) - j dt (N+(t)E AI}M (t) + N+{t)P^ M (t) + A+{t)Q(t)) . 

(2.41) 

Because of that, the total Hamiltonian of the Einstein-Maxwell theory is 

H= [ d 3 x (NH + N S H S + AqQ) + N+(t)E A v M (t) + N+(t)P^ BM (t) + A+(t)Q(t) . (2.42) 

Hence, one is left with the last three terms only when the classical constrainsts 

U = , (2.43a) 

Us = , (2.43b) 

G = (2.43c) 

are satisfied. 

III. BOUNDARY TERMS IN KERR-NEWMAN SPACETIME 

As we saw in the previous section, one must include, in asymptotically flat spacetimes, certain boundary terms in 
order to make the variational principle consistent. Of particular interest are the boundary terms in Kerr- Newman 
spacetime, the most general black hole spacetime. In this section we shall calculate these boundary terms. 

To begin with, we write down the Kerr-Newman line element in Boyer-Lindquist coordinates |27]]: 

j2 A-a 2 sin 2 6> 2 2a sin 2 9(r 2 + a 2 - A) (r 2 + a 2 ) 2 - Aa 2 sin 2 9 . 2 n , , 2 E 2 _ J/l2 ,„ lN 
ds 2 = dt 2 -dtd(j) + sm 2 6d(f) 2 + —dr 2 + Y<d6 2 , (3.1) 

where 

S := r 2 + a 2 cos 2 6 , (3.2a) 
A := r 2 + a 2 + Q 2 - 2Mr . (3.2b) 

In these equations, M is the ADM mass of the hole, Q its charge, and a is the angular momentum per unit mass. To 



calculate the boundary terms we must approximate the line element (3.1) at asymptotic infinity, where r — ► oo, when 
only leading order terms are taken into account: 

2M\ , 2 4Jsin 2 6» ,, 2 . „ „ , (2 A 2M \ , 2 , I/l2 
+2 dido!) + r 2 sin 2 9dd> 2 + 1H Mr 2 + r 2 d9 2 



ds 2 w - 1 oT + r 2 sin 2 6»d<?!) 2 + 1 + dr 2 + r 2 d6» 2 , (3.3) 

\ r J r \ r J 

where J :— Ma is the angular momentum of the hole. In Cartesian coordinates this expression takes the form 

ds 2 w - ( 1 - ^— ] dt 2 - ^{xdy - ydx)dt + ( 1 + - — J (dx 2 + d?/ 2 + dz 2 ) ; (3.4) 
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here r is not the same r as in Eq. (|3.3| ): In Eq. (3.3) r is one of the Boyer-Lindquist coordinates, whereas in Eq. (3.4) 
r is defined to be equal to {x 2 + y + z 2 ) 1 ^ 2 . 

We now proceed to evaluate the boundary terms. When evaluating the boundary terms the first task is to fix 
the coordinate system far away from the black hole. In other words, we must fix the lapse N and the shift N a . In 
this paper we choose a faraway coordinate system which revolves, with respect to the Cartesian coordinates x, y and 
z, with an extremely small angular velocity tu around the z-axis. (We must choose u> to be extremely small since 
otherwise the velocities of faraway observers would exceed the velocity of light.) Because, in flat space, the velocity 
of an observer at the point r — xi + yj + zk revolving with angular velocity Q is 



v = lu x r , 

and because, in Cartesian coordinates, N a represents the a-component of velocity we find that 

iV a = e a bc u: b x c , 



(3.5) 



(3.6) 



where e bc is the Levi-Civita symbol such that £123 = 1. 

What sort of boundary terms do show up with this kind of a choice of the shift? To begin with, we recall from 
Section || that variation of the momentum p ab conjugate to q ab brings along a term 



f) , we must therefore bring along a boundary 



dt / d A x{N a 8p ab )\ b , 



which must be cancelled at infinity. If the shift N a is chosen as in Eq. (|E 
term 



'as 



2e abc / dtu b lim <b x c p an dS, 



(3.7) 



which replaces the boundary term Sq£ 1S of Eq. 

Now, when calculating the boundary t erm Sg% of Eq. (3.7) we should, of course, first perform a coordinate trans- 



formation where the spacetime metric (3A) is replaced by the corresponding expression in revolving coordinates, 
and then proceed to calculate the boundary term by using this expression. However, when the faraway coordinate 
system revolves very slowly, we are interested in terms linear in co only. Taking into account the transformation in 
the expression of the metric wou ld p roduce terms quadrati c in cj, which can be neglected. Hence, we are allowed 
to calculate the boundary term (3.7) by using the metric (3.4). This calculation has been performed in details in 
Appendix and we get, because uj — ujk: 



dtujJ 



(3.8) 



We must still calculate the ADM boundary term ( p.32| ) as well as the electromagnetic boundary term (2.38). The 
ADM boundary term of Kerr-Newman spacetime is, for arbitrary lapse N + at infinity, 



S, 



ADM 

as 



dtN + M 



(3.9) 



To calculate the electromagnetic boundary term we first recall that for Kerr-Newman solution the only non-zero 
components of in Boyer-Lindquist coordinates are 



A, 



Qr 
' E 



. Liar . 2 
Aa, = sin 1 



(3.10a) 
(3.10b) 



Using Eqs. (2.E), (2. IE) and (2.21) one finds that for general lapse and shift one can write p a , the canonical momentum 
conjugate to A a as 



N 4tt 



(3.11) 



This expression, together with Eqs. (3.10), implies that in Boyer-Lindquist coordinates the only surviving component 
of p a is p r which, in the leading order, can be written very far away from the hole as 
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Q 1 0(r- 3 ) . (3.12) 



Anr 2 



Hence, the electromagnetic boundary term (2.3S) is 

' dtA+Q , (3.13) 



as expected. The slow rotation of the asymptotic coordinate system will change the ADM and the electromagnetic 
boundary terms a bit but the resulting corrections will be of the order of 0(u> 2 ) and can therefore be neglected. 



IV. HAMILTONIAN DYNAMICS OF KERR-NEWMAN SPACETIMES 



We shall now proceed to the study of the Hamiltonian dynamics of maximally extended Kerr-Newman spacetimes. 
To begin with, consider a foliation of such spacetimes into space and time. Obviously, we want the spacelike hy- 
persurfaces where the time t — constant to cover as great a portion of spacetime as possible. Maximally extended 
Kerr-Newman spacetimes have a periodic geometrical structure, and we pick up one such period |^7| . We choose the 
spacelike hypersurfaces t = constant such that they begin from the left hand side asymptotic infinity, then go through 
the interior regions of the hole in arbitrary ways, and finally end at the right hand side asymptotic infinity in the 
conformal diagram. However, such spacelike hypersurfaces cannot be pushed beyond the interior horizons where the 
Boyer-Lindquist coordinate 

r = r_ := AI - \J M 2 - Q 2 - a 2 , (4.1) 

since otherwise our hypersurfaces would fail to be spacelike. Hence our study of the Hamiltonian dynamics of Kerr- 
Newman spacetimes must be restricted to include, in addition to the left and the right exterior regions of the Kerr- 
Newman black hole, only such an interior region of the hole that lies between two successive r = r_ hypersurfaces 
in the conformal diagram. Our spacelike hypersurface t = constant therefore begins its life at the past r = r_ 
hypersurface, then goes through the bifurcation point where 

r = r + :=M+ \J M 2 - Q 2 - a 2 , (4.2) 

and finally ends its life at the future r = r_ hypersurface (See Fig. [l]). Bearing this restriction in mind, we shall now 
go into the Hamiltonian dynamics of Kerr-Newman spacetimes. 

The first task is to write the action with appropriate boundary terms. The problem is now that we have two 
asymptotic infinities, and at both of these infinities we have certain boundary terms. When this fact is taken into 
account, we find that the action takes the form 

S = J dt J d 3 x (p ab q ab +P a A a - NH — N S H S - AoG) - J dt [{N+ + N~)M + (A+ - A~)Q + - u~)J] . 

(4.3) 

In this equation, quantities equipped with plus and minus, respectively, are quantities written at the right and the left 
asymptotic infinities. In particular, lo + and lu~ are angular velocities of faraway coordinate systems around z-axis. 
Hence, the total Hamiltonian of Kerr-Newman spacetime is 

H tot = I d 3 x (NH + N S H S + A g) + (N+ + N~)M + (A+ - A^)Q + (lu+ - w~)J . (4.4) 

Now, the problem with our Hamiltonian is that it contains an enormous number of independent degrees of freedom. 
Indeed, our Hamiltonian may be considered as a function of the hypersurface metric q a b at each point x on the spacelike 
hypersurface S, together with the corresponding canonical momenta p ab . However, the ultimate object of interest 
in this paper is canonical quantization of the stationary black hole sector of Einstein-Maxwell theory. Stationary 
black holes, in turn, are characterized by just three independent, classical degrees of freedom, and hence an enormous 
number of degrees of freedom must be truncated. 

For non-rotating black hole spacetimes, the truncation process may be performed in the following manner: One 
first writes the action for asymptotically flat, spherically symmetric Einstein-Maxwell theory. After the Hamiltonian, 
diffcomorphism and Gaussian constraints have been solved, one is left with just four canonical degrees of freedom 
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which can be taken to be the mass M and the electric charge Q of the hole, together with the corresponding canonical 
momenta pu and pq |24|]26| , |28| . A similar truncation could also be performed for rotating black holes: One begins 
with asymptotically flat Einstein-Maxwell theory with appropriate symmetries, solves the classical constraints, and 
is finally left with just six physical, canonical degrees of freedom which may be taken to be the mass M, the electric 
charge Q, and the angular momentum J of the Kerr-Newman black hole, together with the corresponding canonical 
momenta pm, Pq and pj. 

An important feature of the process explained above, in which the phase space becomes reduced in such a way that 
only the physical degrees of freedom are left, is that the resulting Hamiltonian, the so called reduced Hamiltonian, 
involves the boundary terms only. In particular, the reduced Hamiltonian of Kerr-Newmann spacetimes is now 

H rcd = (7V + + N - )M + ( A + _ A -j Q + ( W + _ ^-j j (4 5) 

As a matter of fact, the reduced Hamiltonian may be used as the real, physical Hamiltonian of the system. This was 
proved by Regge and Teitelboim Q . They found that if one solves the classical constraints and then substitutes the 
solutions to the reduced Hamiltonian, then the correct equations of motion for the canonical variables are obtained. 
More precisely, they showed the following: One assumes that the variables q a b and p ab can be separated by a one to 
one, time independent, functionally differentiable canonical transformation in two sets (tp a ,n a ) and (i/j a ,tta) in such 
a way that: 

a) The reduced Hamiltonian depends only on tp a and the n a . 

b) When the 7r a are prescribed as functions p a of x which satisfy 

Pa = 0, (4.6) 

then the constraints Ti = and TL S = can be solved to express the <p a as functionals 

<p a = n^A] (4.7) 



of the remaining canonical variables. 



The functional derivatives of /" with respect to ip A and tta are assumed to exist. If the above conditions are true 
then Hamilton's equation for the reduced Hamiltonian 



H [ip ;tta] = (boundary terms) 



(4.8) 



<p"=f ,TV a =p a 



together with Eqs. (4.6) and (4.7) are equivalent to Einstein's equations in the particular frame defined by ix a = p a . 

The proof of this result is easy: The Poisson brackets are invariant under canonical transformation and the Hamil- 
tonian is unchanged in value if the canonical transformation is independent of time. Hence 



8H 



5tta(x) 



On the other hand, 

H[tp a ;ir a ,il) A ]TTA] 



(boundary terms) 



Differentiating Eq. ( 4.10| ) with respect to tta gives 



d 3 y 



SH 



Sf a (y) 



SH 



5tt a {x) 



However, by Eq. (4.6) 



Tr a (x) 



6H 



8ip a (x) 



= , 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



and therefore 
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6H 5H lcd 

SH „. =/ .„„. " ' «" 3) 

In other words, H lcd generates the correct equation of motion for ip A . In a completely analogous way one shows that 
the correct equation of motion is also generated for it a- Although we have here considered pure gravity only, it is 
clear that our analysis could be easily generalized to include electromagnetic fields as well. 

The real problem is now: Are the assumptions of the previous theorem valid for Kerr-Newman spacetimes? In 
other words, is it possible to divide the phase space of an Einstein-Maxwell theory with appropriate symmetries in 
two parts in a manner explained above? For spherically symmetric, asymptotically flat Einstein-Maxwell theory this 
can be done and has been done in Refs. |^,^|. For theories having the Kerr-Newman solution as their unique solution 
to the classical constraints this has not been done. However, there is not an obvious reason why this could not be 
done, and we state the following hypothesis: 

For an appropriately symmetric, asymptotically flat Einstein-Maxwell theory having the Kerr-Newman solution 
as its unique solution to the Hamiltonian, diffeomorphism and Gaussian constraints there exists a one to one, 
time independent, differentiable canonical transformation which divides the phase space {q a b,P ab , A a ,p a ) into 
two sets (M, Q, J, Pm, Pq, P.j) and (ip A , Pa) in such a way that: 

a) The reduced Hamiltonian depends only on M, Q, J, Pm, Pq and Pj . 

b) When the M, Q and J are prescribed as functions m, q, and l which satisfy 

m = q = i = , (4-14) 

then the constraints can be solved to express the Pm, Pq and Pj as functionally differentiable functionals 
of ip A and Pa- 

We have been unable to find an exact proof of this hypothesis for Kerr-Newmann black hole spacetimes and indeed, 
this is the weak point of our model. However, there are no obvious reasons why it would not be true. In what follows, 
we shall consider our hypothesis as true and see where it takes us. 



The first consequence of our hypothesis is that H lcd of Eq. (kL5|) may be used as the real, physical Hamiltonian of 
our theory, with M, Q and J as the coordinates of the configuration space. For that reason we shall drop "red" from 
our Hamiltonian, and denote it simply by H. 

Our Hamiltonian now implies the following canonical equations of motion: 

OH 

M=- — = 0, (4.15a) 

opm 

' = , (4.15b) 

j = tt- = , (4.15c) 



dp 



j 

P M= -dM=~ {N+ + N ~ h (4 - 1M) 

PQ=-qq=-(K-A q ) , (4.15e) 
Ft H 

pj=-^- = -(u+-u,-) , (4.15f) 

Where pm, Pq and pj, respectively, are canonical momenta conjugate to M, Q and J. As expected, M, Q and J are 
constants in time. The time derivative of pm depends on the choice of the lapse function at both asymptotic infinities 
of our spacelike hypersurface, pq on the difference between electric potentials at infinities, and pj on the difference 
between the angular velocities of faraway coordinate systems. 

The quantities iV^, A and lo^ determine the gauge of our theory. For physical reasons, it is sensible to work in a 
specific gauge where 

N+= 1 , (4.16a) 
N-= , (4.16b) 
w±=0, (4.16c) 
A$=Q. (4.16d) 
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In this gauge the Hamiltonian takes a particularly simple form in terms of the canonical coordinates: 

H = M . (4.17) 

The physical sense of this kind of gauge fixing lies in the fact that we consider Kerr-Newman spacetimes from the 
point of view of a certain specific observer: Our observer is at rest at the right hand side asymptotic infinity, and 
his time coordinate is the asymptotic Minkowski time, the proper time of such an observer. We have "frozen" the 
time evolution at the left infinity, which is sensible because our observer can make observations from just one infinity. 
For such an observer, the classical Hamiltonian of the Kerr-Newman spacetime is just M, the ADM mass of the 
Kerr-Newman black hole. 

Now, the problem with the phase space coordinates M, Q, J, pm, Pq and pj is that they describe the static aspects 
of Kerr-Newman spacetimes only. However, there is dynamics in Kerr-Newman spacetimes in the sense that between 
the event horizons there is a region in which it is impossible to find a timelike Killing vector field. Our next task is 
to find canonical variables describing the dynamical properties of Kerr-Newman black holes in a natural manner. 

When choosing the phase space coordinates, we refer to the properties of our observer: Our observer lies at rest 
very far away from the hole and he is an inertial observer. For such an observer, the Kerr-Newman spacetime appears 
as stationary, and all the relevant dynamics of the Kerr-Newman spacetime is, in a certain sense, confined inside 
the event horizon of the hole. These properties prompt us to choose the phase space coordinates in such a manner 
that when the classical equations of motion are satisfied, all the dynamics is, in a certain sense, confined inside the 
event horizon r = r + of the hole. Moreover, as we shall see in a moment, the choice of the phase space coordinates 
describing the dynamics of spacetime is related to the choice of slicing of spacetime into space and time. We choose 
a slicing where the proper time of an observer in a free fall through the bifurcation surface coincides with the proper 
time of our faraway observer at rest. On grounds of the principle of equivalence one may view these types of slicings 
to be in a preferred position in relating the physical properties of the black hole interior to the physics observed by 
our faraway observer. 



A. Hamiltonian dynamics with charge and angular momentum as external parameters 

To make things simple, consider J and Q first as mere external parameters of the theory. In that case our phase 
space is just two-dimensional being spanned by the phase space coordinates M and pu- In this two-dimensional phase 
space we now perform the following transformation from the " old" phase space coordinates M and pu to the " new" 
phase space coordinates R and pr: 



\p M \= \/2MR -R 2 -Q 2 -a 2 + M sin -1 



M - R 



-7rM 

2 



p R = sgn(p M )V2Mi?-i? 2 -Q 2 -a 2 , 
and we have imposed by hand a restriction 

— ttM < pm < ttM . 



(4.18a) 
(4.18b) 

(4.19) 



With the restriction (4. IE) the transformation (4.18) is well defined and one to one. It follows from Eq. ( 4.18b| ) that 



M = —(p 2 R + R 2 + Q 2 +a 2 ) 



(4.20) 



If one substitutes this expression for M into Eq. (4.18a), one gets pu in terms of R and pr. One finds t hat the 
fundamental Poisson brackets between M and pu are preserved invariant, and hence the transformation (4.15) is 

canonical. 

Equations (4.17) and ( 4.20| ) imply that the classical Hamiltonian takes, in terms of the variables R and pr, the 
form 



h = Jr + r2 + q 



2 ' a 2 ) 



(4.21) 



The geometrical interpretation of the variable R is extremely interesting. We first write the Hamiltonian equation of 
motion for R: 
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and it follows from Eq. (4.21) that when the classical equations of motion for M and pu are satisfied, then the 
equation of motion for R is 



R 



2M 



1 



Q 2 



R 2 



(4.23) 



Now, one can see from the Kerr- Newman metric (3.1) that for an observer falling freely through the bifurcation surface 
at the equatorial plane 9 = n/2 such that 9 = <f> = 0, the proper time elapsed when r goes from r to r + dr is dr such 
that 



dr' 



Q 2 - 2Mr 



-dr 2 



and therefore the equation of motion of our observer is 

r l = 1 

r 



(4.24) 



(4.25) 



where the dot means proper time derivative. As one can see, Eqs. ( 4.23| ) and (4.25) are identical. Hence, we may 
interpret R as the radius of the wormhole throat of the Kerr-Newman black hole, from the point of view of an observer 
in a free f all at the equatorial plane such that <fi — through the bifurcation two-sphere. Moreover, one can see from 
Eq. ( 4.23 ) that R is confined to be, classically, within the region [r_, r + ]. In other words, our variable R "lives" only 
within the inner and the outer horizons of the Kerr-Newman black hole, and this is precisely the region in which 
it is impossible to find a time coordinate such that spacetime with respect to that time coordinate would be static. 
Hence, both of the requirements we posed for our phase space coordinates are satisfied: Dynamics is confined inside 
the apparent horizon and the time coordinate on the wormhole throat is the proper time of a freely falling observer. 
Wi th th e interpretation explained above, the restriction (4. IE) becomes understandable. One can see from 



Eq. ( 4.15d ) that when the lapse functions at asymptotic infinities are chosen as in Eqs. ( 4.16 ), the canonical 
momentum pm conj ugate to M is — t + constant, where t is the time coordinate of our asymptotic observer. Now, 
the transformation (4.18) involves an identification of the time coordinate t with the proper time of a freely falling 
observer on the throat. However, as it was noted at the beginning of this section, it is impossible to push the spacelike 
hypersurfaces t = constant beoynd the r = r_ hypersurfaces in the conformal diagram. The proper time a freely 
falling observer needs to fall from the past r = r_ h ypersurface to the future r = r_ hypersurface through the 



bifurcation surface is, as it can be seen from Eq. (4.25), 



At = 2 



'dr' 



y/2Mr' - r' 2 - Q 2 - a 2 



= 2irM , 



(4.26) 



and hence the restriction ( 4. IE ) is needed. As one can see from Eq. ( 4.18a ), \pm\ = when R — r + and \pm\ = kM 
when R — r_ . We have chosen pm to be positive when the hypersurface t = constant lies between the past r = r_ 
hypersurface and the bifurcation surface, and negative when that hypersurface lies between the bifurcation point and 
the future r = r_ hypersurface. 

Concerning the classical Hamiltonian theory with J and Q as mere external parameters the only thing one still 
needs to check is whether there exists such foliations of the Kerr-Newman spacetime where the Minkowski time at 
asymptotic infinity and the proper time of a freely falling observer at the throat through the bifurcation surface really 
are the one and the same time coordinate. It is easy to see that time coordinates determining this sort of foliations 
do exist. A concrete example is constructed in Appendix^. It should be noted, however, that all foliations in which 
the proper time on the throat and asymptotic Minkowski time are identified are incomplete since such foliations, in 
addition to failing to cover the regions outside the past and the future r = r_ hypersurfaces also fail to cover the 
whole exterior regions of the hole. More precisely, these foliations are valid only when —ttM < t < nM (see Fig. |h. 



B. Hamiltonian dynamics with charge and angular momentum as dynamical variables 



The next task is to complete the classical Hamiltonian (4.21) such that Q and a are replaced by functions of 
appropriate phase space variables describing the dynamics of Kerr-Newman spacetimes in a natural manner. To this 
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end, we must find, for constant M, a canonical transformation from the phase space coordinates (Q,pq) and (J,pj) 
to some new phase space coordinates which we shall denote by u and v, and their canonical momenta p u and p v . 

We shall perform such a transformation in two steps. At the first stage we replace Q and a by canonical momenta 
conjugate to yet some unknown coordinates W\ and W2 of the configuration space: 



Pwi • — Q 



Pu 



and the classical Hamiltonian of Eq. ( 4.21 ) takes the form 

1 



H = 



2R 



R* 



(4.27a) 
(4.27b) 



(4.28) 



The next task is to find w\ and wi . One expects that w\ and u>2 are r elated in one way or another to the momenta pq 
and pj conjugate to Q and J, respectively. Because we see from Eq. ( 4.15e ) that pq determines the electromagnetic 
gauge and from Eq. (4.151) that pj determines the angular velocity of faraway coordinate systems we first write 
the classical Hamiltonian in a general electromagnetic gauge when faraway coordinate systems rotate with arbitrary 
angular velocities: 



H 



1 

2R 



(pr+P 2 Wi +p1 2 +R 2 ) +(A+ -A )p Wl +M(u+ -lu ) Pn 



(4.29) 



which follows from Eq. ( |4.5| ). Using Eqs. ( 4.15c ) and ( 4.151 ) and the fact that M is a constant when the classical 
equations of motion are satisfied, we get the Hamiltonian equations of motion for w± and u>2 



w%:- 



w 2 - 



dH 


Pwi 


dp Wl 


~ ~R 


dH 


Pv)2 


dpw 2 


~ ~R 



PQ 



(4.30a) 
(4.30b) 



An e xpressi on for pp a nd pj in terms of R, pr, wi, u>2, p Wl and p W2 can be gained by integrating both sides of 
Eqs. ( 4.30a ) and ( 4.30h| ) along the classical trajectory in phase space: 



PJ- 



1 

M 



Pw2 j D 

— -an — W2 , 
RR 



where we have substituted 



<2M 

R = -sgn(p M )\/ — - 1 - 



Pi x 



Pw 2 



R 



(4.31a) 
(4.31b) 

(4.32) 



This substitution involves choosing pg = pj = 0. When the electric potentials are assumed to vanish at infinities, 
and the asymptotic coordinate systems are assumed to be non-rotating, this kind of choice can be made. With an 
appropriate choice of the integration constant we get 



PQ-= sga(p M )Pw! 
PJ := sgn(p M " 



sm 



R 2 



. \J (p 2 r + Pli + P' 



2 

wo 



2R 



P\+P 2 W1 



R 2 







K 




| FW2 


sin- 1 | 



AR 2 (p 2 wl +p 2 W2 
Pr + Pi, 




w 1 , 



{Pr+P 2 Wi 



p 2 W2 +R 2 Y -±R 2 {p 2 Wi +pI, 



(4- 
(4- 



where we have made the substitution 



M 



2ft (PR+Pl, 



R 2 



(4.34) 



Equations (4.18b), (4.27) and (4.33) now constitute a transformation from the phase space coordinates M, pm, Q, 
Pq, J andpj to the phase space coordinates R,pr, w\, p Wl , 11)2, andpu, 2 . One can easily show that this transformation 
is well defined and canonical. Moreover, the transformation is one to one provided that we impose the restrictions 
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Mp.J —W2 
P™2 



< 7T 



< 7T 



(4.35a) 
(4.35b) 



These restrictions are related to the fact that we are considering spacetime between two successive r — r_ hyper- 
surfaces. Since both pq and pj vanish when the electric potentials are assumed to vanish at asymptotic infinities 
and the asymptotic coordinate systems are assumed to be non-rotating, we find that classically w\ and w 2 have the 
following properties: At the past r = r_ hypersurface w% — —Qir + pq and w 2 = — an + Mpj, at the bifurcation 
surface w% = pq and w 2 = Mpj, and at the future r = r_ hypersurface w\ — Qn + pq, and W2 = an + Mpj. In 
other words, the classical domains of wi and w 2 are bounded by the fact that the t = constant hypersurfaces cannot 
be pushed beyond the r = r_ hypersurfaces. 

As the last step we perform a canonical transformation from the variables w\, p Wl , u>2 and p W2 to the variables u, 
p u , v and p v ^. We define 



11 . — Pwi '< 



Pu-= Pu 



Pv- 



■ p W2 sin 

Pw 2 C0 S 



Pw± 
Wl 
Pw\ 
W 2 

Pw 2 

W2 

Pw2 



(4.36a) 
(4.36b) 
(4.36c) 
(4.36d) 



This transformation is well defined, canonical and, with the restriction ( 4.35| ), one to one as well. We find that 

(4.37a) 
(4.37b) 



P 2 ■■ 



Pi 



Pw2 ' Pv 



In other words, we may identify p\ + u 2 as the square Q 2 of the electric charge Q, and p 2 + v 2 as the square of a 2 of 
the angular momentum per unit mass of the hole. Because of that, the classical Hamiltonian of Kerr-Newman black 
holes finally takes a very simple form 



Pn 



Pi 



R z 



(4.38) 



V. QUANTUM THEORY OF KERR-NEWMAN BLACK HOLES 



After completing the classical Hamiltonian theory of stationary spacetime containing a Kerr-Newman black hole, we 
are now prepared to consider the canonical quantum theory of such spacetimes. In what follows, we shall concentrate 
on a specific class of canonical quantum theories. More precisely, we define the Hilbert space to be the space 



L 2 (I 



, R s dRdudv) with the inner product 



(*i|*2> := 



dRR s 



du I dv^* 1 (R,u,v)^ 2 {R,u,v) . 



(5.1) 



Through the substitutions pn - 
by the corresponding symmetric operator 



' an 



,Pu 



H := -\K 



J du 



dR \ dR 



and p v 



1 d 2 



we replace the classical Hamiltonian of Eq. (4.3S 



2R8W 2 ~ 2R8V 2 



-R 



u v 
2R + 2R 



(5.2) 



u and ii should not be confused with light cone coordinates or anything like that! 
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This operator may be viewed as the Hamiltonian operator of Kerr-Newman black holes. Its eigenvalues are eigenvalues 
of the ADM energy E of such a hole, from the point of view of a faraway observer at rest. The eigenvalue equation 
in question takes the form: 



~2 R M 



R s 



d_ 

dR 



1 d 2 
2R0V 2 



J__cP_ 

2Rdv i 



1 u 2 v 
2 R+ 2R + 2R 



*(i?, u, v) = E^(r, u, v) 



(5.3) 



This equation is the main result of this paper. In a certain sense, it can be considered as a sort of a "time- independent 
Schrodinger equation of all the possible black holes" , and *&(R , u, v) as the wave function of black holes. Concentrating 
on the quantum theories where s = 1, we find that Eq. (|5.3|) takes a particularly simple and beautiful form: 



1 

2R 



3 2 
dR 2 



du 2 



dv 2 



R l 



*(i?, u, v) = E^(R, u, v) 



If we write 



*(i?, u, v) = 4(R)ip. x (u)ipiv) 



we find that Eq. (5J3) can be separated to eigenvalue equations for M, Q 2 and a 2 

1 



1 , d f , , d 

2 dR \ dR 



:R 



2R 2R 
d?_ 

du 2 

£_ 

dv 2 



i{R) = MiHR) , 
+ u 2 J ipt(u) = Q 2 ipi(u) 
+ v 2 ) f 2 (v) = a 2 (p 2 {v) . 



(5.4) 

(5.5) 

(5.6a) 
(5.6b) 
(5.6c) 



Consider now Eq. ( |5. 6a ), the eigenvalue equation for the ADM mass M of the hole, in more details. It can be 
written as: 



R- s — (r 8 - 1 — 
dR \ dR 



(5.7) 



As one can see, the function 



Q 2 a 2 

li + R +R - 2M 

is negative when r_ < R < r + and positive (or zero) elsewhere. Semiclassically, one may therefore expect oscillating 
behavior from the wave function when r_ < R < r+ and exponential behavior elsewhere. Hence, our system is 
somewhat analogous to a particle in a potential well such that R is confined, classically, between the outer and the 
inner horizons of the black hole. What happens semiclassically is that the wave packet corresponding to the variable R 
is reflected from the inner horizon. As a result, we get, when the hole is in a stationary state, a standing wave between 
the outer and the inner horizons. Thus the classical incompleteness, associated with the fact that our foliation is valid 
only when —nM < t < ttM , is removed by quantum mechanics: In a stationary state there are no propagating wave 
packets between the horizons and our quantum theory is therefore valid at any moment of time. 

When a = Q = 0, we have a Schwarzschild black hole, and the inner horizon is replaced by the black hole singularity: 
The wave packets are no more reflected from the inner horizon but from the singularity. Again, we have a standing 
wave in a stationary state and the quantum theory is valid at any moment of time, but the wave lies between 
Sch warz schild horizon and the singularity. As such there is an interesting resemblance between the properties of 
Eq. ( |5.3|) and those of the Schrodinger equation of a hydrogen atom: When the hydrogen atom is in an s-state where 
the orbital angular momentum of the electron orbiting the proton vanishes, the electron should, classically, collide 
with the proton in a very short time. Quantum mechanically, however, the wave packet representing the electron 
is reflected from the proton, and finally the electron is represented by a standing wave, which makes the quantum 
theory of the hydrogen atom valid at any moment of time. In a Schwarzschild black hole, the proton is replaced by 
the black hole singularity, and the distance of an electron from the proton is replaced by the throat radius R of the 
hole. Nevertheless, the solution provided by quantum theory to the problems of the classical one is similar for both 
black holes and hydrogen atoms. 
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We shall now enter the detailed analysis of the eigenvalue equation (5.6a' 



) . To begin with, we see that if we denote 



and define 



= R 3 ' 2 , 



2s- 1 
6 

7- 2s 



s > 2 



s < 2 



(5.8a) 
(5.8b) 



(5.9a) 
(5.9b) 



then Eq. ( [5.6a ) takes the form 



ax^ x z 



9 



,2/3 



j£x) = Mxtx) 



(5.10) 



This equation has been analyzed in details in Ref. [^9|. The only difference between Eq. ( 5. 10] ) and Eq. (3.18) in 



Ref. jl9) is that Q 2 has been replaced by Q 2 + a 2 . Hence one just replaces Q 2 by Q 2 + a 2 in the results obtained for 
Eq. (3.18) in Ref. @. 

As in Ref. jl4|, one can show that the spectrum of M is discrete, bounded below, and can be made positive. From 
the physical point of view, the semiboundedness and positivity (in some cases) of the spectrum are very satisfying 
results: The semiboundedness of the spectrum implies that one cannot extract an infinite amount of energy from 
the system, whereas the positivity of the spectrum is in agreement with the well-known positive-energy theorems of 
general relativity, which state, roughly speaking, that the ADM energy of spacetime is always positive or zero when 
Einstein's equations are satisfied 



However, one can prove even more than that. One can show that the eigenvalue equation (5.10) implies that when 
p > 3/2, the eigenvalues of the quantity 



M 2 - 



are strictly positive, and when 1/2 < p < 3/2, the eigenvalues of the quantity M 2 — Q 2 — a 2 can again be made 
positive by means of an appropriate choice of the boundary conditions of the wave function yrfx) at the point x = 
or, more precisely, by means of an appropriate choice of a self-adjoint extension. Moreover, the WKB analysis of 



Eq. (|5.10|) yields the result that when Q 2 + a 2 > 1, and M 2 - Q 2 - a 2 > 1 such that r > 1, the WKB eigenvalues 



M n have a property 



Mi 



2n + 1 + 0(1) , 



(5.11) 



wher e n i s an integer and 0(1) denotes a term that vanishes asymptoticall y for large n. A numerical analysis of 
Eq. ( 5.1 0|) yields the result that, up to the term 1 on the right hand side, Eq. ( 5.11 ) gives fairly accurate results even 



when yQ^ + a 2 and n are relatively small (i.e., of order 10). In other words, it seems that the eigenvalues of the 
quantity \J M 2 — Q 2 — a 2 are of the form \[2~n in the semiclassical limit. 



Now, how should we understand these results? The positivity of the spectrum of the quantity M — Q 



a 



has 



an interesting consequence regarding Hawking radiation: If one thinks of Hawking radiation as an outcome of a chain 
of transitions from higher- to lower-energy eigenstates, the positivity of the spectrum of M 2 — Q 2 — a 2 implies that a 
non-extreme Kerr-Newman black hole can never become, through Hawking radiation, an extreme black hole with zero 
temperature, a result that is in agreement with both the third law of thermodynamics and the qualitative difference 
between extreme and non-extreme black holes. One may consider this as a strong argument in favor of our choice of 
the phase space coordinates describing the d ynam ics of Kerr-Newman spacetimes. 

Before considering the implications of Eq. ( 5.11 ), let us calculate the spectra of Q and a from Eqs. ( |5.6b| ) and ( 5.6c ). 
As one can see, both of these equations are, essentially, time -independent Schrodinger equations of a one-dimensional 
linear harmonic oscillator. When the solutions to Eq. ( 5.6b ) are chosen to be harmonic oscillator eigenfunctions, one 
finds that the eigenvalues of Q 2 are 



2k + 1 



(5.12) 



or, in SI units, 
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Ql = {2k + 1) 



(5.13) 



where k = 0, 1, 2, In this equation, e is the elementary charge and 

e 2 1 



Aire^hc 137 



(5.14) 



is the fine structure constant. In other words, Eq. (5^3), the "Shrodinger equation of black holes", implies that the 
electric charge of black holes has a discrete spectrum. 

One may have very mixed feelings on the physical validity of the charge spectrum in Eq. (5.13) 



particles at least, the electric charge Q itself, instead of its square Q is an integer. 



For elementary 
Because of that it might 



appear at the first sight that the charge spectrum we have just obtained contradicts all the possible observations and 
expectations, and should therefore be rejected on physical grounds. 

Such a conclusion, however, would be much too rapid. Firstly, elementary particles are certainly not b lack holes 
because for them \Q\ ^> M. Secondly, a dimensional investigation reveals us that the charge spectrum ( 5.12 ) is exactly 
what one expects for black holes: when one writes the electric charge in terms of the natural constants eo, h and c, 
one finds that the natural unit of electric charge is the so called " Planck charge" 



Qpi := \/ Aneohc 



(5.15) 



One observes that the square Q Pl of the Planck charge Qp\, instead of the Planck charge Qp\ itself, is proportional 
to h. Now, for bounded systems, the observed quantities usually tend to be quantized in such a manner that when 
we write that quantity in terms of the natural constants relevant to the system under consideration, then H must be 
multiplied by an integer in the spectrum. In a hydrogen atom, for instance, the relevant natural constants are eo, fr, 
e and the mass m e of the electron. From these quantities one may construct a natural unit of energy in a hydrogen 
atom: 

4 

m e e 



(^e ) 2 h 2 ' 

and one expects the energy to be quantized such that the energy eigenvalues are of the form 

m „4 



E,, 



"7 



(iTre ) 2 h 2 n 2 ' 



(5.16) 



where 7 is some pure number and n is an integer. Indeed, if we take 7 = 1/2, we get exactly the correct energy 
spectrum for a hydrogen atom. Now, for black holes the o nly natural constants we are allowed to use are, in SI units, 
h, c, G and eo- Hence, the Planck charge Qpi of Eq. ( 5.15 ) is the natural unit of charge for black holes, and therefore 
one expects that t he sq uare of the electric charge, instead of the charge itself, must be an integer. In other words, the 
charge spectrum (5.13) is exactly what one expects for black holes. 

Let us now turn our attention to Eq. (5.6c) which gives the spectrum of a 2 . As for the electric charge, we find that 
the possible eigenvalues of a 2 are 



= 21 + 1 



or, in SI units, 



where I = 0, 1, 2, . . .. Again, 
is multiplied by an integer. Putting Eqs 
eigenvalues of the black hole are 



or, in SI units 



9 /~ , 1 KG 
a 2 = (21 + 1) — 

c 



(5.17) 



(5.18) 



one observes that the quantity under consideration is quantized in such a way that K 
(5.12) and (5.17) together we find that, semiclassically, the mass 



(5.11 



(5.19) 



/ 2mMpi , 



(5.20) 
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where 



m := n + 1 + k = 0, 1,2, . 



(5.21) 



and 



Mi 



PI ■— \ T7 



(5.22) 



is the Planck mass. 

The spectra of the quantities M, Q and a now have interesting consequences regarding the area spectrum of black 
holes. As it is well known, the area of the outer horizon of a Kerr-Newman black hole is 



whereas the area of the inner horizon is 



A =An(r 2 _+a 2 ) 



Using Eqs. ( 5.11 ) and ( 5.13 ) we observe that the semiclassical eigenvalues of the quantity 

A tot := A+ + A- , 

which we shall call, for the sake of convenience, the total area of a black hole, are of the form 

16?r(2n + 21 + k) 



A tot 



or, in SI units, 



where 



A tot 



167r(2n + 2Z + fc^p! , 



hG 



£pi := \l — — 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



(5.28) 



is the Planck length. In other words, we have obtained a result which is closely related, although not quite identical to 
the proposal suggested by Bekenstein in 1974 and since then revived by several authors: According to that proposal 
the spectrum of the outer horizon of the black hole is of the form f|-|25| 



A+ = in l 



pi i 



(5.29) 



where n is integer and 7 is a pure number of order one. Hence, we have obtained a result w hich states that the total 
area of the hole, with 7 = 167r, instead of the area of its inner horizon, is quantized as in Eq. (5.2E). In the last section 
of this paper we shall consider in more details the possibility that it is perhaps the total area, and not the area of the 
outer horizon, which should be an integer in Planck units. 

As the final check of our quantum theory of black holes, let us calculate the angular momentum spectrum of black 
holes. We observe from Eqs. ( [5.18 ), (5.20) and (5.21) that the possible eigenvalues of the angular momentum J = Ma 
of the hole are, semiclassically, of the form 



J n , hk - ±2y/l(l + n + k)h 



(5.30) 



For uncharged black hole where k = we therefore find, in the limit of extremality where I 3> n, that the angular 
momentum eigenvalues are of the form 



rrijh 



(5.31) 



where rrij = 0,±2,±4, 



As one can see, the angular momentum spectrum of black holes, as predicted by our theory is, at least in the limit 
of extremality, exactly what one might expect. Even the fact that the angular momentum J is an even number, is in 
harmony with our expectations: When the black hole performs a transition from one angular momentum eigenstate to 
another, a graviton is emitted or absorbed. Because the spin of a graviton is two, one might expect that the angular 
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momentum of the black hole could change only by an even number. For instance, one may show, quite rigorously 
that when a system consisting of two mass points revolving around their common center of mass emits or absorbs a 
graviton, the angular momentum quantum number of t he sy stem can change only by an even number p9fl . Because 



of that, the angular momentum spectrum given by Eq. (5.31) for extremal black holes may be used as a very strong 



argument in favor of the physical validity of our quantum-mechanical model of black holes. 

Unfortunately o ur m odel also appears to contain a very serious problem regarding the angular momentum spectrum: 
According to Eq. ( |5.30 ) the angular momentum of a black hole is not in general an integer times the Planck constant 



H. Should we be worried because of this result? 

The answer to this question is: Not necessarily. The usual rules for the quantum mechanics of angular momentum 
follow from the symmetries of flat spacetime, and spacetime containing a Kerr-Newman black hole is certainly not 
flat. In curved spacetime the angular momentum eigenvalues of a system do not necessarily have the same properties 
as they would have in flat spacetime. 

To illustrate this fact by a simple example, consider a particle moving in a conelike spacetime geometry (See Fig. 3). 
The z-component L z of the angular momentum eigenvalues may be calculated from the equation 

- iK-g7%4>) = L z %4>) , (5.32) 
from which it follows that the angular momentum eigenfunctions are of the form 

= Ce^ L '* , (5.33) 

where C is a constant. In flat spacetime the period of \P((/>) is 2ir, producing the usual angular momentum spectrum. 
In conelike spacetime geometry, however, the period of ^ is not 2tt but 2tt — e, where e is the deficit angle of the cone 
(See Fig. ^|). In other words, we must have: 

+ 2tt - e) = , (5.34) 
and therefore the angular momentum eigenvalues are of the form 

L z = m z ^ , (5.35) 

where m z = 0, ±1, ±2, .... In other words, the angular momentum of a system is not necessarily an integer times the 
Planck constant in curved spacetime. 

VI. CONCLUDING REMARKS 

In this paper we have considered a particular quantum-mechanical model of Kerr-Newman black holes. The fun- 
damental ideas behind our model were based on the black hole uniqueness theorems. According to these theorems a 
black hole in stationary spacetime is completely characterized by exactly three free variables which may taken to be 
the mass M, the electric charge Q and the angular momentum J of the hole. From these theorems it follows that the 
Kerr-Newman solution, being completely characterized by these three free variables, is the most general stationary 
black hole solution to combined Einstein-Maxwell equations. In our model we considered a Hamiltonian quantum 
theory of stationary black hole spacetimes in such a way that the phase space was spanned by the variables M, J 
and Q, together with the corresponding canonical momenta pm, Pj and pq. The problem with these phase space 
coordinates, however, is that they describe the static aspects of black hole spacetimes only. However, there is dynamics 
in Kerr-Newman spacetimes in the sense that between the horizons there is no timelike Killing vector field, and we 
managed to find new phase space coordinates which describe the dynamical properties of Kerr-Newman spacetime in a 
particularly natural manner. These phase space coordinates were replaced by the corres pon ding quantum mechanical 



operators yielding the symmetric Hamiltonian operator. Our analysis produced Eq. (5.3) which, in a certain very 
restricted sense, may be considered as a sort of "Schrodinger equation of black holes". That equation gives, in the 
context of our model, the mass, electric charge and angular momentum spectra of black holes. 

Eq. ( |5.3] ), which is the main result of this paper, implies that the mass, electric charge and angular momentum 
spectra of black holes are discrete. Moreover, it implies that the mass spectrum is bounded from below and can be 
made positive. By means of a choice of an appropriate self-adjoint extension one may show that the spectrum of the 
quantity 

M 2 -Q 2 -a 2 , 
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where a is the angular momentum per unit mass of the hole, is always positive. Regarding Hawking radiation, this is 
a very important result: It means that a non-extreme black hole can never become an extreme black hole by means 
of the Hawking radiation of black holes. This result is in agreement with the third law of black hole thermodynamics, 
and is therefore a strong argument in favor of the physical validity of our model. 



At the high end of the spectrum, Eq. (5.3) implied that the eigenvalues of the quantities M, Q and a are all 
quantized, in natural units, in a very similar manner: In natural units the eigenvalues of these quantities are all of 
the form \/2n, where n is an integer. Although this kind of a spectrum might appear very odd for an electric ch arge 
spectrum of black holes, it is exactly what one expects on dimensional grounds. In the extremal limit, Eq. ( |5.3| ) 
implied that the angular momentum eigenvalues of black holes are of the form 

rrijh , 

where rrij — 0, ±2, ±4, . . .. 

Of particular interest is the area spectrum of black holes given by Eq. (5.3). Eq. (5.3) implied that the sum of the 
areas of the two horizons of Kerr-Newman black hole is of the form 

where Zpi is the Planck length. Hence, we get a result which is closely related, although not quite identical to, the 
proposal made by Bekenstein in 1974. According to Bekenstein's proposal, it is not the sum of areas of horizon but 
the area of the outer horizon which has an equal spacing in its spectrum. 

Although our result about an equal spacing for the spectrum of the sum of the horizon areas may have certain 
esthetic merits, it also involves some problems. For instance, the fact that the mass eigenvalues are of the form \j2m 
which, together with the fact that Q and a have similar spectra, implied the area spectrum under consideration, also 
implies that the angular frequencies of quanta of Hawking radiation emitted in transitions between nearby states is 

LO ~ — — . (6.1) 

M V ' 

For Schwarzschild black holes this is something one might expect because the Hawking temperature of such a hole 
is @ 

T„ = — !— . (6.2) 

and therefore it follows from Wien's displacement law that the maximum of the thermal spectrum of black hole 
radiation corresponds to the angular frequency 

Wmax OC . (6.3) 

In other words, the angular frequency associated with the discrete spectrum of Hawking radiation as predicted by our 
model, behaves, as a function of M, in the same way as does the angular frequency corresponding to the maximum 
of the thermal spectrum as predicted by Hawking and others. 

Unfortunately, this nice correspondence between Hawking's result and our model breaks down when Q or a are 
different from zero. In that case the Hawking temperature of the black hole is (30| 

r„ = . V«>-Q>-°> , , , (6 . 4) 



2tt 



(m + v/Af 2 - Q 2 - a 2 ) +a? 



and one finds that the maximum of the thermal spectrum corresponds to the angular frequency 

VM 2 - Q 2 - a 2 



M + v/Af 2 - Q 2 - <r 



(6.5) 



In other words, the angular frequency (6.1) predicted by our model corresponds, when the hole is near extremality, 
to a temperature which is much higher than the Hawking temperature. 

However, there may be a possible way out of this problem. In all our investigation we have emphasized the 
importance of the dynamics of the intermediate region between the horizons of the black hole. The dynamics of this 
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intermediate region is, in our model, responsible for the discrete eigenvalues of the mass, electric charge and angular 
momentum of the hole. Now, if we take this point of view to its extreme limits we are tempted to speculate that 
both of the horizons of the hole, acting as the boundaries of the intermediate region, may participate, in one way or 
another, in the radiation process of the black hole. In other words, both of the horizons may radiate. The radiation 
emitted by the inner horizon is probably emitted inside the inner horizon, and is therefore not observed by the external 
observer. Nevertheless, an emission of this radiation is likely to reduce considerably the number of quanta, and hence 
the temperature, of the radiation coming out from the hole: The more the inner horizon radiates, the less quanta are 
left for the outer horizon. 

Let us give up for a moment our resistance to this most charming temptation and have a play with the thought 
that both of the horizons have an important role in black hole radiation. For instance, one might consider one quarter 
of the total area of the hole as a sort of a "total entropy" of the hole 

Stot = ~(A+ + A-) . (6.6) 

Moreover, one might be inclined to define a temperature T corresponding to this entropy (whatever that means) such 
that 



1 dStot ,„ -■. 

T ~dE~ ' (6 - ?) 



and one finds, quite unexpectedly, that 



In other words, we have recovered the Hawking temperature of the Schwarzschild black hole (see Eq. ( J6 . 2[ ) ) . This 
expression is the same for all black holes, and it is inversely proportional to the mass M of the hole. It may well be 
that all this is just meaningless play with symbols, without any physical content, but nevertheless the idea that it 
is the total area, and not the area of the outer horizon, which is of fundamental importance in black hole quantum 
mechanics, appears to possess remarkable internal consistency: If the total area of the hole has equal spacing in its 
spectrum, one expects the temperature of the hole to be inversely proportional to the mass M, and this result is 
recovered if the total entropy of the hole is taken to be one quarter of not the area of the outer horizon but of the 
total area of the hole. We shall investigate these ideas in more details in forthcoming papers. 

To conclude, our quantum-mechanical model of Kerr-Newman black holes appears to involve several physically 
sensible properties but also some problems. For instance, the claim that Kerr-Newman spacetime and our phase 
space variables satisfy the assumptions of Reggc's and Tcitclboim's theorem has been left unproved. The proper 
analysis of the Hamiltonian dynamics of Kerr-Newman spacetimes along the lines shown by Kuchaf for Schwarzscild 
spacetime should therefore be performed ||(J . 

Another problem is, whether the quantum mechanics of black holes can be described with a sufficient accuracy by 
means of a model having just three independent degrees of freedom. In other words, are the mass, electric charge 
and angular momentum spectra obtained from our model reliable? When answering to this question one can just say 
that at least the spectra are such as one might expect on semiclassical and dimensional grounds. As to the problems 
related to the statistical origin of black hole entropy and things like that our model says nothing. Another, more 
esoteric, reason why our model may probably contain some hints of truth is its simplicity and certain naturality. Such 
things, however, are merely matters of taste and should not be trusted too much. 
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APPENDIX A: A BOUNDARY TERM 



In this Appendix the boundary term So£ of Section III is calculated explicitly. 



That boundary term is expressed in Cartesian coordinates in Eq. (3.7). Because the coordinate system revolves 
around the z-axis with angular velocity u>, we have to 1 = = lu 2 and u> 3 — u>, and so the boundary term can be 
written in the form 
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S&g = -2 J dt ui I (xp 2 s dS s - yp\dS s ) 



(Al) 



By comparing the line element of Eq. (3.4) with the ADM line element 

ds 2 = — (N 2 — N a N a ) dt 2 + 2N a dx a dt + q ab dx a dx b , (A2) 

where N is the lapse function, N a is a component of the shift vector (a = 1,2,3) and q a b is a spacelike component of 
the metric tensor associated with the hypersurface, we find that the only non-zero spacelike components of the metric 
tensor are 



2M 

qu = H = 922 = <733 

r 



and for the components of the shift vector N a we have 



2Jx 



N 3 = 



where 



and 



N a := q ab N b 



r = \/ x 2 + y 2 + z 2 



Comparing Eqs. (3.4) and (A2) we obtain the lapse function N: 

/ 2M 4J 2 / 2M\ _1 " " ~ 
N=\l + -g- 1+ (x 2 +y 2 ). 



(A3) 

(A4) 

(A5) 
(A6) 

(A7) 



To evaluate covariant derivatives of N a we must calculate Christoffcl symbols using the non-zero components of q a b 
expressed in Eq. (A3). The non-zero Christoffcl symbols needed in the calculations are 



Mx 



iu = r? 2 = r? 3 = --=-! + — 



2M 



v 1 - r 1 - 

1 22 — 1 ■« — 



My 



, , 2M 

rii r l! = — (i + — 

2M 
r 

r 2 -r 1 -r 3 + — 

1 22 — 1 12 — 1 23 ~ „3 1 ' „ 



33 



Mx 

t-3 



2M 



r 1 

1 13 



- r 2 - 

— 1 23 — 



By using the formula 



Kab = ^ (-«0i + ^a|6 + ^6|a) 



for the exterior curvature tensor K a b, we get 



K n = 



2M 



2Jxy 
Nr 5 W + 2M 



- 3 



K l2 
K 13 
K 23 



K 22 
K33 
K21 

K31 

K 32 



2M 



2Jxy 

' Nr 5 \r + 2M 



= , 



Nr 5 



(- 2 - v 2 ) ( 



2M 
r + 2M 



- 3 



2M 



Jyz 

Nr 5 \r + 2M 



Jxz 



2M 



Nr 5 \r + 2M 



(A8a) 
(A8b) 
(A8c) 
(A8d) 
(A8e) 

(A9) 

(AlOa) 

(AlOb) 
(AlOc) 
(AlOd) 

(AlOe) 

(AlOf) 
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To evaluate p a b, the canonical momentum conjugate to q a b, from the expression 

Pab ■= TF~ ( K ab ~ QabK) , 
±07T 

we must first calculate K and ^/q, and we get 

K := q 11 !^ + q 22 K 22 + q 33 K 33 = , 



2Af\ 3 



r / 



When the results from Eqs. (A12) and (A13) are substituted to Eq. (All), we have 



1 + 2M) 2 2 Jxy ( 2M 



Pu 



P22 



16tt Nr J \r + 2M 
{l + ^) i 2Jxy( 2M 



- 3 



16tt iVr 5 Vr + 2M 



3 , 



P33 = , 



Pl2 = P21 



Pl3 = P31 



P23 = P32 



(1 + ^) 7 J 

16?r Nr 5 



2Af 
r + 2M 



16tt TVr 5 V^ + 2M 
fl + ^-Y Jxz ( 2M 



16tt Nr 5 \r + 2M 



(AH) 

(A12) 
(A13) 

(Al4a) 

(A14b) 
(Al4c) 

(A14d) 
(Al4e) 
(A14f) 



The boundary term Sq£ can be expressed by using the components calculated above and the components of the unit 
normal on the surface: 



'as 



2 / dt u 



xq (p2W +P22n +p 2 3n ) dS - yq {pun +p\2n + pi 3 n J 



dS 



1 

8^ 



dt loJ (h 1 



2M 



V 



2M 



r ) \r+2M 



Nr 5 



(x 2 — y 2 ) (xn 1 — yri 2 ) + 



2xy (yn 1 — xn 2 ) + z(x 2 + y 2 ) 



2\ 3 



dS . 



(A15) 



This integral is easy to evaluate in spherical coordinates. We first consider a 2-dimensional spherical surface with 
radius r. The relations between the spherical coordinates r, 9 an (j) and the Cartesian coordinates x, y and z are 



x = r cos <p sin 0, y = rsmasm<p, z = rcost 
The components of the unit normal n a ,(a = 1,2,3) on the surface are 



n = n x = cos</>sin#, n 2 — n y — sin 9 sin (j), n 3 — n z — cost 



and the area element is 

dS = r 2 sin 

In these coordinates the boundary term takes the form 



(A16) 
(A17) 
(A18) 



'as 



1 

8^ 



dt oj J 1 



2M 



V 



2M 



r J \r+2M 



- 3 



1 

N 



cos 



sin 2 <h) 2 sin 5 



8=0 <ji=0 



4 cos 2 6 sin 2 (A sin 5 9 + cos 2 sin 3 I 



(A19) 
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As r approaches infinity N goes to 1 and so the denominator can be approximated as 1. Integration gives then: 

and so the boundary term at infinity, where r — > oo, is 

Sgg = - j dtuJ . (A21) 



APPENDIX B: A NOVIKOV-TYPE SLICING OF KERR-NEWMAN SPACETIME 

In this appendix we construct in details a slicing of Kerr-Newman spacetime in which the time coordinate of a 
freely falling observer through the bifurcation surface and the flat Minkowski time of a faraway observer at rest at 
the right hand side asymptotic infinity are identified. In a certain sense, one may view these observers and their time 
coordinates as physically equivalent. In Refs. H] and similar identifications are performed and they are based 
on the Novikov coordinate system (see, for instance, Ref. where the time coordinate of a given point is given 

by the proper time r of a freely falling observer in the Schwarzchild or Reisner-Nordstrom spacetime through that 
point, and the radial coordinate R* in the Novikov coordinate system is related to the point r where the freely falling 
observer has begun his journey. 



Since the ^-coordinate in the classical Hamiltonian ( 4.21 ) can be geometrically interpreted as the radius of a 



wormhole throat at the equatorial plane 8 — tt/2 in the Kerr-Newman black hol e, w e begin the construction of 
the slicing with desired properties by considering the Kerr-Newman line element ( p.l[ ) written in Boyer-Lindquist 
coordinates at the equatorial plane: 

, 2 A -a 2 , 2 2a(r 2 + a 2 - A) , , , (r 2 + a 2 ) 2 - Aa 2 ,, r 2 , , 

ds = la - dt 3 + 3 d ^> + -£ dr i ( B1 ) 



where 



A = r 2 + a 2 - 2Mr + Q 2 . (B2) 
When the Boyer-Lindquist coordinates x M (/i = 0, 1, 2, 3) satisfy the constraint 



yJ-g^xW = 1 , (B3) 

where g^ u gives the components of the metric tensor of the Kerr-Newman spacetime, the Lagrangian of a particle in 
the Kerr-Newman spacetime is, in general, 

£kn = -^g^x^x 11 , (B4) 



At the equatorial plane Eqs. (Bl) and (B4) give for the Lagrangian an expression 



lA-a 2 , a(r 2 + a 2 - A) . • 1 (r 2 + fl 2 ) 2 - > 2 .„ lr 2 
= + " 2^ - 2 A r ' (B5) 

where the dot denotes proper time derivative. The canonical momenta conjugate to t, r and <fi are: 

9Lkn A -a 2 . a(r 2 +a 2 ~ A) 
Pt = »■ = n — t H 4> , (B6a) 



dt 

Pr = -g— = --^r , (B6b) 



8L KN (r 2 + a 2 ) 2 - Aa 2 ■ q(r 2 +a 2 -A) ; 

■ = ~ </>H t 2 1. (Bbc) 



It is easy to see that 
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Pi 



dt 

d(j> 



, 
. 



Therefore the momenta pt and p^ are constants of motion. Let us denote these constants as 

R* 

Pt 



P4> =■ £ , 

where we have introduced new real valued parameters R* and 
The Hamiltonian 



Hkn = p r r + p t t + P4,(j> - Lkn 
of a particle in Kerr-Newman spacetime can be shown to coincide with the Lagrangian Lkn: 



Therefore, when the constraint (B3) is satisfied, we get 



On the other hand, when we use the parameters introduced in Eqs. (B8), we get from Eq. (BE): 



2Hi 



R* 



KN 



Vl + R* 



A + , 



and from Eqs. (B6t) and (B6c) we obtain tp and t in terms of I and R*: 



1 r A + 2r A a 2 -a A (A-a A 



R* 



a(r 2 + a 2 - A) 



A \ r 2 

1 / a(r 2 + a 2 - A) 



R* 



Vl + R* 2 
A -a 2 



VI + R* 



Now, we choose £ = 0, and because of this particular choice we get from Eq. (B12) 

1 (r A + 2r 2 a 2 -a 2 (A-a 2 ) 



2H } 



R* 



R* 



KN 



A VT+r^ a 



Vl + R* 



As we set r = 0, Eq. (B14) yields us a quartic equation for 



Ar 2 = 



R* 



1 + R* 



[r 4 + 2r 2 a 2 -a 2 (A-a 2 



(B7a) 
(B7b) 

(B8a) 
(B8b) 

(B9) 
(BIO) 

(Bll) 
(B12) 



(B13a) 
(B13b) 

(B14) 

(B15) 



From this equation one can calculate the r-coordinate r mm of the point from which an observer in a free fall begi ns hi s 
journey, in terms of R* which will henceforth be used as a radial coordinate of Kerr-Newman spacetime. Eq. (B14) 
implies an implicit expression r(r, R*) for the "old" radial coordinate r in terms of the "new" time coordinate r and 
the "new" radial coordinate R*: 



±Vl + R* 2 [ 

J r 



,a*(ii*) ^R* 2 ( r > 2 + a 2 ){2Mr' - Q 2 ) - r' 2 A 



dr' 



(B16) 



In this equation the signs + and — , respectively, correspond to the past and the future of the line where the time 
coordinate t = i n the conformal diagram. To obtain an explicit expre ssion r(r,R*) f or r one should first solve the 
quartic equation ( B15 ), and then perform the integration in Eq. ( |B16| ). Solving Eq. ( B15 ), however, would yield a 



2G 



tremendously complicated expression for r max , and we shall not write it down here. However, it is easy to see that 
ther e are always at least two positive roots r = r max = r max (i?*). This can be seen by plott ing the both sides of 
Eq. ( B15| ) and varying R* . Moreover, o ne fin ds that if one puts r = r max = r + then Eq. ( |B15 ) implies R* — 0, and 
vice versa: if one sets R* = 0, then Eq. ( B15 ) is solved by r = r+. Hence, we have found that for every R* > there 
is an observer in a free fall such that this observer is at rest at the time t = r = with respect to the "old" radial 
coordinate r. When R* = our observer begins his journey at the bifurcation surface and his world line is a straight 
vertical line in the conformal diagram. 

Can we extend this coordinate transformation to the right hand side asymptotic infinity? Yes we can, since we may 
choose the coordinate R* such that the solution r = r 



> U 



choice is made, one can show starting from Eq. (B15), that for large r 



is the largest of the roots of Eq. ( B15| ). When this 



- 1 



Q 2 - 2Mr n 



+ <?(r-L) • 



(B17) 



Hence R* goes to infinity as r max goes to infinity and vice versa. Moreover, the time coordinate r of an observer at 
the asymptotic infinity coincides with the proper time r of a freely falling observer at the wormhole throat. 

Another matter to investigate still is that do the observers rotate or not with respect to the Boyer-Lindquist 
coordinates? We wrote our Hamiltonian from the point of view of an asymptotic non-rotating observer, and we 
assumed a foliation in which the time coordinate at the throat is a proper time of a non-rotating observer in a free 
fall. To show that in our foliation both of the observers are non-rotating we must show that <fi — > as r — > r + and 
r — > oo. The latter case is straightforward, since in the expression 



1 a(r 2 



R* 



VT+R* 



(B18) 



given by Eq. (B13b) when I = 0, the factor 



R* 



— — ^ approaches to one and the factor in front of it approaches to zero. 
The first case where r — > r + is a bit tricky, since we do not know the explicit relation of r and R* at the bifurcation 



point. We have solved the tricky part by expanding Eq. (B15) in terms of r near the bifurcation point. If we take only 
the zeroth and the first order terms, we find that the point r = r max where r — is related to R* by an expression 



r + I 2 



R* 2 



l+R* 



2Mr% + o 2 ri™r + AMa 2 r 2 ^r 
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2r: 



2Mr% + 2Ma 2 r +J ^ - 2a 2 q 2 j^ 



(B19) 



l+R* 



which gives that r max = r + as R* — > 0, as it should. Now, when Eq. (B19) is substituted into Eq. ( B1S| ) and R* — > 0, 
one gets the result 







(B20) 



In other words, we have managed to construct a foliation of Kerr-Newman spacetime with desired properties at the 
equatorial plane: At the asymptotic infinity the time coordinate is the proper time of a freely falling, non-rotating 
observer at rest, and at the wormhole throat that of a similar non-rotating observer in a radial free fall through the 
bifurcation surface. 

It is even possible to show that our construction gives the Novikov coordinate system in the Schwarzschild spacetime 
when one sets q = a = in Eq. (B15). This result is given by Eqs. (B15) and (B16). We get an analogous coordinate 
system for the Reissner-Nordstrom spacetime when only a = 0. It can be shown that then the relation between r max 
and R* is 



M 



M 2 - q 2 (l + R* 2 ) (l + R* 
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FIG. 1. The conformal diagram of Kerr-Newman spacetimes. Our spacelike hyp ersurf aces t = constant begin their life at the 
past r = r_ hypersurface, then go through the bifurcation point, and finally end their life at the future r — r_ hypersurface. 




FIG. 3. A particle moving in a conelike spacetime geometry. When the cone is stretched on a plane, the deficit angle 
e appears. As a result of the appearance of this deficit angle, the periodic boundary condition for the angular momentum 
eigenfunction ip{<j)) is tft(<p + 2n ~ e) = 
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